Astronomy & Astrophysics manuscript no. 5370 


© ESO 2008 


February 5, 2008 





A closure model with plumes 
II. Application to the stochastic excitation of solar p modes 

o : 

K. Belkacem 1 , R. Samadi 1 , M. J. Goupil 1 , F. Kupka 2 , and F. Baudin 3 

Or 

m : 
o\ : 



Observatoire de Paris, LESIA, CNRS UMR 8109, 92195 Meudon, France 
Max-Planck-Institute for Astrophysics, Karl-Schwarzschild Str. 1, 85741 Garching, Germany 
Institut d'Astrophysique Spatiale, CNRS/Universite Paris XI UMR 8617,91405 Orsay Cedex, France 



Received 06 April 2006 / Accepted 30 June 2006 

ABSTRACT 



CO 

> 
o 

, Context. Amplitudes of stellar p modes result from a balance between excitation and damping processes taking place in the uppermost part of 

■ convective zones in solar-type stars and can therefore be used as a seismic diagnostic for the physical properties of these external layers. 
Aims. Our goal is to improve the theoretical modelling of stochastic excitation of p modes by turbulent convection. 

■ Methods. With the help of the closure model with plume (CMP) developed in a companion paper, we refine the theoretical description of the 
' excitation by the turbulent Reynolds stress term. The CMP is generalized for two-point correlation products so as to apply it to the formalism 

*^ I developed by Samadi & Goupil (2001). The excitation source terms are then computed with this improvement, and a comparison with solar 
Q i data from the GOLF instrument is performed. 

Results. The present model provides a significant improvement when comparing absolute values of theoretical ampplitudes with observational 
data. It gives rise to a frequency dependence of the power supplied to solar p modes, which agrees with GOLF observations. It is shown that 
the asymmetry of the turbulent convection zone (up and downflows) plays a major role in the excitation processes. Despite an increase in the 
Reynolds stress term contribution due to our improved description, an additional source of excitation, identified as the entropy source term, is 
still necessary for reproducing the observational data. 
j^J ' Conclusions. Theoretical excitation rates in the frequency range v 6 [2.5 mHz, 4 mHz] now are in agreement with the observational data from 
k>( | the GOLF instrument. However, at lower frequencies, it exhibits small discrepancies at the maximum level of a few per cent. Improvements are 
. likely to come from a better physical description of the excitation by entropy fluctuations in the superadiabatic zone. 

. p.: 
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1. Introduction 

Amplitudes of solar-like oscillations result from a balance 
between excitation and damping. Excitation is attributed to 
turbulent motions that excite the p modes. In the uppermost 
part of the convection zone, entropy fluctuations and eddy 
motions drive oscillations. In this region, convection becomes 
inefficient and there is an increase of in the eddy velocities and 
entropy fluctuations. Solar-like oscillations are mainly excited 
in such a region, thus a theoretical model of the excitation 
processes is a powerful tool in understanding the properti es of 
the co nvective zones of solar-type stars. IGoldreich & Keelevl 
l l 1977b have proposed a model for the excitation process using 
the turbulent Reynolds stress and deducing an estimation of 
the power supplied to the p modes. The underestimation of 
the excitation rates by around a fa ctor 10 3 compared to the 
observed solar values (Osaki 1990) led to alternative formu- 
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lations (Goldreich & K umarl Il990t IGoldreich et ail Il994l) 
Anoth er source of excitation was identified bv lGoldreichetall 
( 1994J): the so-called entropy sou rce term. Its contr i bution 
cannot be neglected, even though IStein & Nordlundl feOOll) 
have shown that excitation from the Reynolds stress remains 
dominant in comparison with the entropy fluctuation source 
term. 

ISamadi & Goupill J200ll) propose a generalized formalism, 
taking the Reynolds and entropic fluctuation source terms 
into account. This approach allows investigation of the effects 
of several models of turbulence ( Samad i et alJ l2003blal) by 
expressing the source terms as functions of the turbulent 
kinetic energy spectrum and the time-correlation function. 
A confrontation of t his model with data from the BiSON 
instrument (data from lChaphn et alJll998) led to the conclu- 
sion that the theoretical predictions were in g ood agreement 
with the observations JSamadi et alJ l2003ah Nevertheless, 
observational data from the GOLF instrument and a study 
of the BiSON data indicate that some discrepancies remain 
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between the theoretical computation and observational data. In 
ISamadi & Goupill J200ll) (see also Samadi et al., 2005), one 
of the main assumptions is the quasi-normal approximation 
(QNA), which is useful for correlation functions of the turbu- 
lent Reynolds stress and the entropy fluctuation source terms 
ISamadi & GoupifcOOlh . 

The uppermost part of the convection zone being a turbulent 
convective system composed of two flows, the probability 
distribution function of the fluctuations of the vertical velocity 
and temperature does not follow a Gaussian law (Lesieur 
Il997h . Thus, the use of the QNA, which is exact for a normal 
distribution, becomes a doubtful approximation. 
In a companion paper ( Bel kacem et alJl2006i hereafter Paper 
I), we propose another approach in order to build a closure 
model that expresses fourth-order correlation functions in 
terms of the second-order ones. This alternative approach 
consists in considering the convection zone as composed 
of t wo flows (the updraf t s and downdrafts). Starting from 
the Grvanik & Hartmann (2002) approach, we develop a 
generalized two-scale mass-flux model (GTFM) that takes the 
physical properties of each flow into accoun t. Then a theoret- 
ical d escription of the plumes developed by iRieutord & Zahnl 
is used to construct the closure model with plumes 
(CMP). This model is valid for one-point correlation functions 
and in the quasi-adiabatic zone. However, what is needed here 
is a closure model for two point correlation functions. In the 
present paper, we then propose a simple way to obtain this 
closure mode l to use it for calcu lating of the excitation rates 
according to Samadi & Goupill J200 lh . Only the Reynolds 
stress source t erm is corrected, mainly because it is th e 
dominant term dStein & Nordlundl200lUSamadi et al.ll2003al) . 
The entropy fluctuations are considered in the same way as 
explained in Sam adi & Goupill J200ll) (i.e. using the QNA 
approximation). 

The paper is organized as follows: the theoretical model of 
stochastic excitation of p modes is briefly summarized in Sect. 
2. In Sect. 3, the closure model with plume (CMP) is gen- 
eralized for two-poi nt correlation products an d implemented 
into the formalism of Samadi & Goupi| J200ll) . In Sect. 4, the 
calculation of theoretical power is explained. In Sect. 5, GOLF 
observational data are presented together with the derivation of 
observable quantities. A comparison between the theoretical 
power and heights computed as described in Sect. 4 with 
the corresponding observed quantities defined in Sect. 5 is 
performed in Sect. 6. Section 7 is dedicated to discussions and 
conclusions. 



2. A model for stochastic excitation of solar-like p 
modes 

The theoretical mo del of stochastic ex citation considered here 
is basically that o f lSamadi& Goupill fcOOl ) (see also Samadi 
et al., 2005). It takes two sources into account that drive the 
resonant modes of the stellar cavity. The first one is related to 
the Reynolds stress tensor and as such represents a mechanical 
source of excitation. The second one is caused by the advec- 
tion of the turbulent fluctuations of entropy by the turbulent 



motions (the so-called "entropy sourc e term") and thus rep- 
resents a thermal sourc e of excitation (Goldr eich et aljfl994l: 
IStein &Nordlundl200ll) . 

The power fed into each mode, P, is given by (see e.g. 
ISamadietalfeOOll) : 



P = ^ =2/7£ = 77<|A| 2 )/^ 



(1) 



where () denotes the ensemble average, (|A| 2 ) the mean square 
amplitude, r\ the damping rate, and E the energy that is defined 

as 



E = 1 <|A| 2 > Ioj 2 



(2) 



where / is t he mode inertia and cap is the oscillation eigenfre- 
quency (see lSamadi & Gonpil200ll for det ails). 
The m ean square amplitude, as explained in lSamadi & Goupill 
(2001), is 



(3) 



where C\ and C 2 . are the turbulent Reynolds stress and entropy 
contributions, respectively. Their expressions for radial modes 
are given by 

d 3 x p 2 f r dre~^ ld 3 r(w 2 wl) (4) 

<J —CO U 

d\ ogr dTe-^ T 



d\ ((wsdx (ws,) 2 > (5) 



where w is the vertical component of the velocity, s, the tur- 
bulent entropy fluctuation and f r (£ r ,m) = f^J , where is 
the radial component of the eigenfunction, and g r a function 
that i nvolves the first and s econd derivatives of (see Eq. 
(9) of lSamadi etail J2003bl) ). Quantities labelled with 1 and 
2 denote two spatial and temporal positions, hence (wjW^ 
and {(ws t )\ {ws t )2 ) correspond to two-point fourth-order cor- 
relation products. These correlation products are usually ap- 
proximated by express ions involving second-o rder products 
only (closure model). In lSamadi & GoupilfeOOll) . the simplest 
approximation was used i.e the quasi-normal hypothesis. We 
study here consequences of using a closure model closer to re- 
ality (i.e the CMP from Paper I). Both are recalled in the next 
section. 

3. Closure models 

3.1. The quasi-normal approximation 

The QNA JLesieurl Il997l Chap VII-2) is adopted in 
ISamadi & Goupill (|200 lb as a convenient means to decompose 
the fourth-order velocity correlations in terms of a product of 
second-order vertical velocity correlations, that is, one uses 



{w\w\)q NA = 2(w 1 w 2 ) 2 + {w\){w\) 
((ws,) A (ws,) 2 ) QNA = (wiw 2 ) (s n sa) 



(6) 



where s t is considered as a passive scalar. 

This approximation (Eq. lHJl) remains strictly valid for normally 
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distributed fluctuating quantities with zero mean. As sho wn by 
iRraichnanl Jl957h in the context of turbulent flows and Stein 
( 1967) in the solar context, the cumulant (the deviation from 
the QNA) can be large and therefore not negligible. The CMP 
presented in Paper I was shown to be a significant improvement 
on the QNA for the one-point correlation products. However, 
we need two-point correlation products here (see Eqs. © and 
l|5}). A generalization of the CMP for two-point correlation 
products is therefore developed in Sect. l3.2l below. 
The second-order correlation products in Eq. are expressed 
in the Fourier domain (k, a>) where k and u> are the wavenum- 
ber and the frequency a ssociated with a turbulent element (see 
ISamadi & Goupill200lL for details). 

3.2. The closure model with plumes 

The closure model with plumes (see Paper I) has been estab- 
lished only for one-point correlation products. Here we gener- 
alize the CMP to two-point correlation products. We start in 
Fig. [Dby comparing the correlation product < w\w\ > calcu- 
lated directly from 3D numerical simulations obtained from the 
Stein & Nordlund code (see Sect.© with those calculated using 
Eq. l|6]of the QNA with second-order correlation products taken 
from the 3D simulation. The question is whether the modelling 
of the k dependency on the two-point correlation function by 
the QNA can be used. For the sake of simplicity, we assume 
that the QNA can be used for the u> dependency. 

The correlation products {w\w^} in Fig. [2 are normalized 
so as to compare only the k dependency of these quantities. 
In the quasi-adiabatic region, the line width at half-maximum 
of the QNA and the numerical product are roughly the same. 
Discrepancies at high values of AX (the correlation length) are 
expected to have a negligible influence on the correlation prod- 
uct. Hence, we assume that the modelling of the k dependency 
on the two-point correlation product by the QNA is valid due 
to a small difference between the line width at half-maximum. 
Hence it is legitimate to use the (k, to) dependency given by the 
QNA. One then needs only to correct the value of the correla- 
tion product at (k — 0, u> = 0) (which corresponds to the one 
point correlation function) with the CMP (see Paper I) for the 
turbulent Reynolds stress term contribution. We use the inter- 
polation formula of iGrvanik & Hartmannl J2002I) for the FOM 
of the velocity (Paper I, Eq. (13)) 

< w\w\ > C mp= (1 + ^Sl) < w\w\ > QNA , (7) 

with {w\w\)qna given by Eq. (jfji the skewness S K is calculated 
from the CMP (see Paper I for details). 

In Fig. |3 calculations using Eqs. (|6jl and are com- 
pared to the direct numerical correlation product. The above 
generalization of the CMP to two-point correlation products 
provides a good approximation mainly in the quasi-adiabatic 
region where the CMP is the more accurate one (see Paper I). 
The k dependence is approximatively modelled by the QNA 
(Fig. ^ except for large correlation lengths (AX > 0.2 Mm), 
but these contribute only negligibly to < w\w\ >. However, in 
the superadiabatic zone, the generalization of the CMP and the 
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Fig. 1. Fourth-order correlation function calculated in the 
quasi-adiabatic zone directly from the 3D numerical simulation 
(solid line) and using the QNA approximation (Eq. (|6j; dashed 
lines). The fourth-order moments are presented as a function of 
the correlation length (AX), and the two curves are normalized 
so as to emphasise only their k dependency. 



QNA both fail to describe the two-point correlation function. 
In that zone, the temperature gradient is varying quickly, which 
is not the case in the CMP. In the plume model (Paper I) the 
temperature gradient appears only through a polytropic law, 
and for sake of simplicity we assume an isentropic atmosphere. 
In addition, for mode lling the FOM < w 4 > , the i nterpolated 
formula derived by IGrvanik & H artmannl d2002» (Paper I, 
Eq. (13)) is not valid in the superadiabatic zone. Thus, in this 
zone the treatment of Eqs. (|6jl and will introduce an energy 
excess injected into high-frequency p modes. 



4. Calculation of the theoretical p mode excitation 
rates 

The rate (P) at which energy is injected per unit time into a 
mode is calculated according to the set of Eqs. © - © when 
the QNA is used and Eqs. ©-Q using the CMP (see Sect.l4.1>. 
The calculation thus requires the knowledge of four different 
types of quantities: 

1) Quantities that are related to the oscillation modes: the 
eigenfunctions (£ r ) and associated eigenfrequencies (u>o). 

2) Quantities that are related to the spatial and time-averaged 
properties of the medium: the density po, the vertical veloc- 
ity w, the entropy S, and a s = dPo/ds. 

3) Quantities that contain information about spatial and tem- 
poral correlations of the convective fluctuations: E(k), 
E s (k), and^(w). 

4) Quantities that take anisotropies into account: a and <1>. The 
value of a is the mean horizontal fractional area of the up- 
drafts (see Paper I), whereas O measures the anisotropy of 
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Fig. 2. Fourth-order correlation function calculated in the su- 
peradiabatic zone (at the top) and in the quasi-adiabatic zone 
(at the bottom) directly from the 3D numerical simulation (dot- 
ted line), using the QNA approximation (Eq. |6j dashed lines) 
and using the CMP (Eq.0 solid line). 



turbu lence and is denned ac cording to lGoughldl977l) (see 
also lSamadi & Goupill200ll for details) as: 



(D = 



< w > 



(8) 



where u 2 = w 2 + u 2 h and m/, is the horizontal velocity. 



Both a and <D are necessary to describe the flow because 
a measures the geometric anisotropy between up and down- 
flows while corresponds to the measure of the velocity 
anisotropies. However, these two quantities are linked because 
of mass conservation. An explicit relation can be easily derived 
between them using the formalism developed in Paper I to ob- 
tain 



a(l - a)5w + a < w > u +(l - a) < w >d 
a(l - a)5w 2 + a < u 2 > u +(\ - a) < u 2 >d 



where the " refers to the velocities of only one flow (updraft 
or downdraft) and 6w is defined as in Paper I. For consistency 
reason, a and O are provided by the 3D numerical simulation. 

4.1. The solar case 

Calculations of the eigenfrequencies and eigenfu nctions (in 
point 1) above) are performed as in Sam adi et alJ d2Q03bh on 
the basis of a ID solar model built according to Gough's (1977) 
non-local formulation of the mixing-length theory (GMLT 
hereafter). 

The spatial and time-averaged quantities in point 2) are ob- 
tained from a 3D simulation of the solar surface. The 3D sim- 
ulations used in this wo rk were built w ith Stein & Nordlund's 
3D num erical code ("see lStein & Nordlundll998llSamadi et alJ 
20033)- Two simulations with different spatial mesh grids are 
considered, namely 253x253x163 and 125x125x82, in order 
to verify that the results are not sensitive to the spatial mesh 
resolution. 

Finally, for the quantities in point 3) the total kinetic 
energy contained in the turbulent kinetic spectrum, E(k), its 
depth dependence, and its ^-dependence are obtained directly 
from a 3D simulation of the u ppermost part of the solar 
convective zone. It was found in ISamadi et alJ J2003al) from 
3D simulations that a Gaussian - usually used for modelling 
Xk - is inadequate: a Lorentzian fits the frequency dependence 
of Xk best. Hence, we adopt a Lorentzian here for^. 



4.2. Calculation of the power injected into the solar 
p modes with the CMP 

We use the generalized CMP for two-point correlation func- 
tions presented in Sect. 13.21 (Eq. Q) to model the Reynolds- 
stress source term. By replacing Eq. with Eq. J7J in Eq. @, 
the calculation of C 2 , (as in Samadi &Goupill2 001) yields: 



fi 15 Jo 



dm{\ + -S w )p 



~r) Jo 



dk 



E 2 (k ) 
k 4 



(10) 



(9) 



Equation (II Oi shows that the CMP causes an increase in the 
power injected into p modes in comparison with calculation 
using only the QNA. On the other hand, the entropy source 
term, C 2 , is still co mputed using the QNA closure model (see 
ISamadi & Goupill200ll for details). 



5. Observational data and inferring observed 
excitation rates 

The observational data set selected here for comparison with 
theoretical calculations was o btained with the G OLF instru- 
ment, onboard SOHO. GOLF JGabriel et aill 19971) is a spec- 
trometer measuring velocities of the photosphere integrated 
over the whole solar disc. Its location on the space platform 
yields a very good signal-to-noise ratio and also continuous 
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Fig. 3. Rate P at which acoustic energy is injected into the 
solar radial modes. Only the Reynolds stress c ontribution is 
compu ted. Cross dots represent P computed from lBaudin et alj 
J2005h solar seismic data from the GOLF instrument (see 
Sect. |5jl. The associated error bars take into account uncertain- 
ties both from the line width (F v ) and from the maximum height 
of the mode profile (H). The curves represent theoretical values 
of P computed as explained in Sect.0] dash-dotted lines corre- 
spond to the calculation of P using the QNA closure model, and 
solid lines represent P using the CMP for the Reynolds stress 
term. We present the results in linear (at the top) and logarith- 
mic scale (at the bottom). 



observations (the actual duty cycle reaches almost 100%). This 
latter characteristic greatly improves the signal to noise ratio in 
the Fourier spectrum. 

However, GOLF suffers from some technical problems, which 
restricts the measurements to one wing of the Na Dj line in- 
stead of both wings. This results in a more difficult absolute 
calibration of the measured velocity and thus a possible bias 
(which does not exceed 20% in terms of the acoustic rate of ex- 
citation). Characteristics of the data set used here are described 



10 7 
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v [mHz] 

Fig. 4. Mode height H calculated as explained in Sect.]5\using 
only the Reynolds stress contribution. The solid (resp. dash- 
dotted) line represents H calculated with the CMP (resp. QNA) 
closure model, and cross-dots represent GOLF data with as- 
sociated error bars. Error bars associated with the curves are 
due to mode line widths that are taken from observations (see 
Eq. O). 



in lBaudin etaflfcOOSl) . 

These observations correspond to two periods when GOLF was 
observing in the same instrumental configuration (blue wing 
of the Na line) with a duration of 805 and 668 days, starting 
on April 11, 1996 and November 11, 2002, respectively. The 
level of solar activity was different during these two periods, 
but the measured excitation rate shows no dependence on ac- 
tivity, as the increase in width compensat es for the decr ease 
in height of the peaks, a s shown by Chan lin etalJ |2000) or 
I.Timenez-Reves et alJ d2003h . 

The GOLF results were compared to BiSON observations and 
are compatible with them over a wide frequency range. A dis- 
crepancy appears at high frequency (v > 3.2 mHz). As the 
height and width of peaks in the Fourier spectrum are af- 
fected by the presence of noise and gaps in the data (see 
Chaplin et al 2003), GOLF was chosen for the comparison 
model/observations. We consider only the € = 1 modes for 
which their properties (line-width, amplitu de) are more accu- 
rately determined than the I — modes (see lBaudin et alJ2 005. 
for details). 

In order to compare theoretical results and observational 
data, the mode excitation rates are inferred from the observa- 
tions according to the relation 



P obs (w ) = 27rr v Mv2( Wo ) 



(11) 



where M = I/^{h) is the mode mass, h the height above the 
photosphere where oscillations are measured, F v = rj/n the 
mode linewidth at half maximum (in Hz), and v] the mean 
square of the mode surface velocity. The last is derived from 
the observations according to 

v 2 = nHr v C ohs (12) 
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where H is the maximum height of the mode profile in the 
power spectrum and C b s the multiplicative constant factor tha t 
depends on the observation technique (see Bau din et all20 05). 
Equation dl 2l i supposes that the mode line profiles are symmet- 
ric, but it is well kno wn that the mode pr ofile deviates from a 
Lorentzian. However, Baudin et al show that this equa- 

tion is accurate enough for the evaluation of the mean square 
of the mode velocity, Eq. dl 2i . On the other hand, the mode 
asymmetry is taken into account when determining mode line 
widths from observational data. 

The mode mass is very sensitive to altitude at high fre- 
quency (see Fig. 1 of Baudin et al 2005), so the layer (h) where 
the mode mass is evaluated must be properly estimated to de- 
rive correct values of the excitation rates. Indeed, solar seismic 
observations in Doppler velocity are usually measured from a 
given spectral line. The layer where oscillations are measured 
then depends on the height where the line is formed. The GOLF 
instrument uses the Na I D 1 and D 2 lines whose height of for- 
mation is estimated at the height h * 340 km (see Baudin et al. 
120051) . 

As an alternati ve to compar i ng the oretical results and ob- 
servational data, Chap lin et alJ (J2005) propose to derive the 
maximum height of the mode profile (H) from the theoretical 
excitation rates and the observed mode line width according to 
the relation: 



27r 2 Mr2C obs ' 
where C a b s = 2.59 for ( — 1 modes. 

Representation of the excitation rates themselves (Eq. (II 1» 
emphasises disagreement at high frequencies, whereas dis- 
agreement at low frequency is more apparent with a representa- 
tion using the profile height (Eq. H3\ ). Note that in the case of 
the observable height, only the slopes are the meaningful quan- 
tities, as the amplitude magnitude depends on the phase of the 
solar cycle when the observations were recorded. 

As the maximum height H strongly depends on the obser- 
vation technique, one cannot compare values of H coming from 
two different instruments. In Fig.|6] we therefore plot the prod- 
uct //C bs, a quantity that is less dependent on the observational 
data (but still through M). Note that for ease of notation, HC b s 
is noted H in the following. 

It is important to stress that the mode height (H) calculated 
from the theoretical excitation rates (Eq. dl It ) depends on the 
observations through the line width F v . This is why in Figs. [3] 
and[6]error bars appear in the theoretical results. In any case, the 
observational data can be characterised by at least three main 
features that the theoretical calculations (see above) must re- 
produce: 

1 . the frequency dependence from low to medium frequencies 
(v < 3 mHz). 

2. the maximum of amplitude at 3 mHz for H and the slope for 
frequencies between 3 and 4 mHz or a nearly flat maximum 
between v a; 3.8 mHz and 4 mHz for P. 

3. the slope at very high frequencies v > 4 mHz. 



6. Comparison between theoretical and observed 
excitation rates 

6.1. Turbulent Reynolds stress contribution 

Figure [3] compares the observed power P injected into solar 
p modes with the theoretical one computed with only the tur- 
bulent Reynolds stress term assuming either the CMP or the 
QNA closure models. Figure 4 shows the associated heights H 
as computed according to Eq. (13). The comparison shows that 
the closure model has a significant effect on the resulting exci- 
tation rates. Indeed, the CMP induces an increase in the energy 
injected into the mode by about a factor two in comparison with 
the QNA closure model and brings the theoretical excitation 
rates closer to the observational ones. This energy increase is 
not uniform in terms of frequencies, due to the variation in the 
skewness with the depth (z) (see Paper I for details) and to the 
fact that the mean square velocity amplitudes of the turbulent 
elements decrease with depth. Indeed, at the top of the convec- 
tion zone where the highest frequency modes are confined, the 
inefficiency of the convective transport causes an increase in 
the velocities. Thus the effect of the flow anisotropy becomes 
dominant for such high-frequency modes. 
At low frequencies (v < 2.5 mHz), the turbulent Reynolds 
stress contribution reproduces the observed power P (Fig. [3} 
within the observational uncertainties. As best emphasised in 
Fig. 4, it is possible that the theoretical results are slightly over- 
estimated, although this remains within the observational error 
bars. 

At intermediate frequencies 4 > v > 3 mHz), the turbulent 
Reynolds stress term is not sufficient to reproduce the observa- 
tions, so the additional excitation coming from entropy fluctu- 
ations is necessary. 

At high frequencies v > 4 mHz), Observational data seem to 
indicate a decrease in the power, which is not reproduced by 
the theoretical power. 

6.2. Adding the entropy fluctuation contribution 

To proceed further, we add the C| contribution (Eq. l|5}). 
Results for the excitation rate and the maximum height are pre- 
sented in Figs.[5]and[6] respectively. The additional (positive) 
entropy contribution causes an overall increase in the excitation 
rates as shown in Fig. [5] The theoretical modelling now repro- 
duces the maximum of the power supplied to the modes when 
compared with the observational data. For the frequency be- 
haviour of the excitation rate and height, (Figs.[3]and[6j show: 

At low frequency (v e [1.6 mHz; 3 mHz]). We pointed out 
in Sect l6.1l that the contribution from the Reynolds stress term 
can be sufficient for reproducing the GOLF data, perhaps even 
overestimating it. The combination of both Reynolds stress and 
entropy fluctuation is too large compared with the observation, 
and the resulting slope differs from the observational one in 
this frequency domain. Note however, that in Fig.[5]error bars 
represent 1 cr error bars. (Fig.|5|i. 

For intermediate and high frequencies (v e [3 ; 4 ] mHz), 
the Reynolds (CMP) and entropy excitation model reproduces 
the v variation in P. This is confirmed with the H representa- 
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tion (Fig.|5]at the bottom). However from a theoretical point 
of view, the description of the behaviour at high frequencies 
(v > 4 mHz) is more complicated because these p modes are 
mainly excited in the superadiabatic zone, which is difficult 
to model properly. On the observational side, it must be kept 
in mind that even data with a signal-to-noise ratio as good as 
GOLF lead to linewidths difficult to measure at high frequen- 
cies. 

7. Discussions and conclusions 

We use a closure model (CMP, Paper I) that is more realistic 
than the usual QNA approximation to model the correlation 
products in a semi-analytical description of the excitation pro- 
cess of solar p modes. The present excitation model gives the 
theoretical slope of the power at intermediate and high frequen- 
cies (v e [2.5 mHz; 4 mHz]), which agrees with the observed 
data. We also find that including the CMP causes a global in- 
crease in the injected power. This brings the power computed 
with the Reynolds stress contribution alone closer to (although, 
at intermediate frequency, still below) the observations. On the 
other hand, the power obtained by including both the Reynolds 
stress and the entropy fluctuation contributions reproduces the 
observations at the maximum of the excitation rates. The com- 
parison can now be made in linear scale, hence at lower fre- 
quencies there is still a small over-estimation (which amounts 
roughly to a few per cents and the errors bars represent lcr 
error bars). The reason for this overestimation cannot be at- 
tributed to the CMP. Indeed, the Reynolds stress contribution 
was compared to the 3D numerical simulation (see Paper I), 
and the one-point fourth-order moment < w 4 > was found to 
agree with the simulation result. The remaining departure from 
the numerical simulation shows that the CMP actually under- 
estimates the FOM in the quasi-adiabatic region, so correcting 
for this bias would result in an even larger overestimation of 
the power. 

Various sources of discrepancies are likely to exist: the sep- 
aration of scales used in the formalism that consists in as- 
suming that the stratification and the oscillations have char- 
acteristic scale lengths larg er than the edd ies contributing to 
the excitation (see Sa madi & Go unil 2001, for details). The 
physical description of the outer layers in the ID solar model 
can also play an important role directly through the veloc- 
ity and indirectly through the eigenfunctions. In this paper, 
we use Gough's (1977) non-local formulation of the mixing - 
lenght theory which shows an improvement in comparison with 
the local formulations in terms of the maximum of power P 
JSamadi et alJl200r3) by about a few percent. Concerning the 
excitation model itself, some improvements in the modelling 
of Reynolds and entropy contributions that ought to be investi- 
gated are outlined below. 

7.1. Turbulent Reynolds stress tensor contribution 
shortages 

At low frequencies, a possibly small overestimation of the 
Reynolds stress contribution can be attributed to the frequency 
dependent factor (xk, see Eq. 10 in Sect. 14. H . IChaplin et alJ 



(2005) use a Gaussian Xk whereas Sam adi etall J2003bl) use 
a Lorentzian factor. In Fig. |6] we present the calculation as- 
suming a Gaussian and a Lorentzian for xu- As shown there, 
the frequency-dependent factory is likely between these two 
regimes. In the quasi-adiabatic convection zone, plumes are 
well-formed, and the convective system must be treated as 
composed of two flows (see Paper I). Hence, the upflows 
that are less turbulent can be modelled by a white noise 
(Gaussian), but downflows are turbulent creating a departure 
from a Gaussian. We expect this effect to cause a decrease in 
the theoretical power and bring it closer to the observation. A 
rough idea can be obtained by taking this effect into account 
as follows: we split the computation of the power supplied into 
the modes into two parts. Those parts correspond to upflow (x^. 
Gaussian) and to downflows (xk- Lorentzian) . The result indi- 
cates a decrease in the power at low frequency, which brings the 
theoretical power closer to the observation. This is true mainly 
for low-frequency modes, which are less sensitive to the su- 
peradiabatic zone where plumes are formed, because this re- 
gion cannot be modelled by such a simple model. This issue 
needs further investigation. 

7.2. Entropy source contribution shortages 

In the present model, the turbulent entropy fluctuations are as- 
sumed to behave as a passive scalar, in other words, the entropy 
fluctuations are assumed to be advected by the turbulent veloc- 
ity field without dissipation. It means that the entropy field does 
not have any effect on the velocity field. 

This assumption associated with the QNA has the advan- 
tage of simplifying the closure of the fourth-order moments 
involving the entropy fluctuations (see Eq. j3.U ). However the 
biases introduced by this assumption remain to be evaluated. 
If the biases turn out to be large, alternative models must be 
developed. 



7.3. Perspectives 

Finally, we stress that there is an additional dependency, the 
coefficient a, which is the mean fractional area of updraft 
on the horizontal plane (see Eq.9). It is a measure of the 
asymmetry of the flows and a small variation in its value 
plays a major role on the excitation rates. This parameter 
has been fixed here using the results of 3-D simulations. 
The influence of parameter a is very important, as a small 
variation of its value leads to an increase in power P through 
the skewness S w (see Paper I). It is beyond the scope of this 
paper to estimate the true effect of a variation in this parameter 
because its value is linked to the physical properties of the 
flows through, for instance, conservation of the mass flux. 
Hence a consistent approach is to investigate a set of different 
numerical simulations. 

The CMP closure model, indeed, strongly depends on the 
structure of the upper convection zone, which again empha- 
sises that the structure of this region is very important in the 
theoretical prediction of the power injected into the p modes, 
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because the skew introduced by the asymmetry increases with 
the departure of a from the value 0.5. It is then possible to 
obtain physical constraints on the asymmetry of the convection 
zone flows. 

To understand what can affect a is therefore an important 
issue, and in near future it will be necessary to study the 
variation in a with the type of star and from a hydrodynamical 
point of view to determine what the main processes that are 
responsible for this asymmetry. One interesting issue is the 
influ ence of a magne tic fi eld on this paramet er: as described 
bv lWeiss etail J2002h and lVoderet"atl(l2005l) . the effect of a 
strong magnetic field induces a reduction in the typical length 
scale of convection, as well as the structure of the flows (hence 
the value of a). 

The study of the mean fractional area a as a function of 
the magnetic field intensity therefore represents an inter- 
esting perspective for characterising B from the excitation 
rates, at least for stars with an expectedly strong magnetic field. 
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Fig. 6. Mode height H calculated as explained in Sect, fus- 
ing only the Reynolds stress contribution. Solid lines represent 
H calculated with the CMP closure model and dots-line is the 
same except that a Gaussian is used for^. Crosses represent 
GOLF data with associated error bars. Error bars associated 
with the curves are due to mode line widths which are taken 
from observation (see Eq. H3\ ). 



Fig. 5. Top: Rate (P) at which acoustic energy is injected 
into the solar radial modes as a func tion of frequency. C ross 
dots represent P computed from the Bau din et alJ (|2005) so- 
lar seismic data from the GOLF instrument (see Sect. [3}. The 
curves represent theoretical values of P computed as explained 
in Sect. 0] the solid line represents P using both the Reynolds 
stress (using the CMP) and entropy source contributions. The 
dotted line corresponds to the calculation for the Reynolds 
stress term using only the CMP closure. Boffom.Mode height 
(H) calculated as explained in Sect. [5] The solid line repre- 
sents H calculated with the CMP closure model, using the 
Reynolds stress and entropy fluctuation contributions. The dot- 
ted line represents H computed with the CMP closure model, 
using only the Reynolds stress contribution. Cross-dots repre- 
sent GOLF data with the associated error bars. Error bars as- 
sociated with the curves are due to mode line widths that are 
taken from observations (see Eq. ( I13» . Only observations near 
minimum solar activity have been used, and they correspond to 
the second period as explained in Sect. [5] 



